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We study a pair of anharmonic microwave cavities that is connected by an optical fiber. The 
photonic spectral density characterizes the evolution of the coupled cavities after the system has 
been prepared in a Fock or NOON state. We evaluate the photonic spectral density within the 
recursive projection method and find that the anharmonicity has no substantial effect on the spectral 
properties. In all cases the photonic spectral density has a Gaussian envelope for large photon 
numbers. 

PACS numbers: 42.50.Ar, 42.50.Pq, 42.50.Ct 



I. INTRODUCTION 



The experimental preparation of Fock states in a microwave cavity [H, 0| is a crucial step towards a 
systematic study of correlated many-body systems with photonic states. Once a Fock state has been 
prepared, we can change the parameters of the system and study the evolution of the Fock state within 
the Hilbert space that is provided by the eigenstates of the Hamiltonian of the new system. This can be 
realized, for instance, in an experiment with two microwave cavities. The system is prepared such that 
one cavity has N photons, the other cavity has no photon. Then the initial Fock state |iV, 0) is assumed 
to be the groundstate of a Hamiltonian Hq. At the time i = the Hamiltonian is switched to the 
Hamiltonian H . This problem has been studied intensively for atomic systems, using the Hubbard model 
and related models [3|-[8l • In the case of two microwave cavities this switch can be realized by connecting 
the two isolated cavities with an optical fiber. Then Hq is the Hamiltonian of the isolated cavities, and H 
is the Hamiltonian of the connected cavities. This type of system, including atomic degress of freedom, 
was studied recently within a Hartree-Fock approximation '9] . In this paper we shall study the dynamics 
of photons which tunnel between two microwave cavities through an optical fiber (cf. Fig. [l}. The effect 
of atoms inside the cavities is approximated by an anharmonic photonic term. A central question is the 
effect of this anharmonicity on the evolution and whether or not the spectral properties depend on the 
initial conditions. 



II. DYNAMICS OF AN ISOLATED SYSTEM 



We consider a system which is isolated from the environment. In terms of photonic states this can 
be realized by an ideal microwave cavity. Furthermore, we assume that the system lives in an A'^ + 1- 
dimensional Hilbert space. With the initial state j^E'o) we obtain for the time evolution of the state 



or the evolution of the return probability |(^'o|'I't)P with the amplitude 
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-iHt\ 



(1) 



FIG. 1: Two microwave cavities coupled by an optical fiber. The energy levels are indicated as those of anharmonic 
cavities. 
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A Laplace transformation relates the return amplitude with the resolvent through the identity 

(*o|*t) - ^(*o|(z - H)-^\^o)e'"''dz , 

where the contour F encloses all the eigenvalues Ej {j — 0,1,..., A^) of H. With the corresponding 
eigenstates \Ej) the spectral representation of the resolvent provides a rational function: 

N 



3=0 

where P]^(z), QN+iiz) are polynomials in z of order N, iV + 1, respectively. These polynomials can 
be evaluated by the recursive projection method (RPM) [10] . This method is based on a systematic 
expansion of the resolvent (5'o|(z — H)~^\'^o), starting from the initial state |^'o). The method can also 
be understood as a random walk in Hilbert space. However, there is a crucial difference in comparison 
with the conventional random walk because the RPM has a prefered direction: The walk of the RPM 
visits a subspace 'H2j only once and never returns to it. A step from 'H2j to 'H2j+2 is defined by the 
Hamiltonian H in such a way that 'H2j+2 is created by acting H on 112] ■ In terms of A'^ photons, 
distributed over two microwave cavities, these subspaces are given by the basis {|A^ — J: j)) bi ^ j)}) 
and the a step is an exchange of a single photon. Thus, the walk follows a path with increasing numbers 
j. This is the main advantage of the RPM that allows us to calculate the resolvent (vPoK^: — ^f)~^|5'o) 
on a A^ + 1-dimensional Hilbert space exactly. In case of the photons in two microwave cavities the RPM 
is efficient for values of A^ up to the order of w lO"'. 

The expression in Ec^. ^ suggests the introduction of the photonic spectral density (PSD) Pc{E) as 
the imaginary part of the resolvent: 

p,{E) = /m(vI.o|(£; - ze - H)-^\^,) = ej^-M^^^^ ■ (3) 

j=Q ^ ^ ^ j) 

The amplitude of the return probability then reads as the Fourier transform of the PSD 

(*o|*t) = lim / p,{E)e-'^'dE . (4) 

Moreover, the overlap (vI'o|i?j) between the initial state and the eigenstates of the Hamiltonian H are 
directly obtained at the poles of the PSD. In the following the PSD shall be evaluated with the RPM, 
using Fock states |0, A^), |A^, 0) and superpositions of these states. 

III. SINGLE CAVITY 

We begin with the simplest case, a harmonic cavity which has a harmonic oscillator spectrum 

with cavity mode ojq and the ceation (annihilation) operator (a) of a photon. The eigenstates are Fock 
states with A^ photons |A) (A = 0, 1, ...) and energy En (A|i?o|A) = ujqN. 

A. Single atom: Jaynes-Cummings model 

An anharmonicity can be created by adding atoms to the cavity which interact with the photons. After 
introducing a single two-level atom to the cavity we can describe the absorption and emission of photons 
by the atom approximately with the Jaynes-Cummings model [11], whose Hamiltonian reads 

Hjc — LJoo'a + {ujQ + A)c'c — g{a^c + c^a) . 
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A is the detuning between the atomic excitation energy and the photon energy, and (c) is the creation 
(annihilation) operator of the atomic excitation. The eigenvalues of this Hamiltonian are (Tl| 

En± = ^o(^ + 1/2) ± + V(iV + 1) (5) 

with eigenstates 

|7V,±) = i=[±|iV,l) + |iV + l,0)] . 

Thus, the eigenstates are not Fock states but a superposition of two Fock states, one with N photons 
and the atom in the groundstate and one with N — 1 photons and the atom in the excited state. The 
transition amplitude from the initial Fock state {N, 0|\I't) is a simple oscillation because only |7V — 1, ±) 
contributes. The superposition of many Fock states can lead to a more complex behavior, such as a 
collapse and revival dynamics [H, [l^ . Ignoring the atomic degrees of freedom, the effect of the atom on 
the photonic spectrum is the additional square-root term in Eq. 



B. Many atoms: anharmonic cavity 

The presence of many atoms in the cavity is quite complex due the presence of the atomic degrees 
of freedom. Therefore, we apply a simplified picture by ignoring the atomic degrees of freedom and 
describing their effect on the photons as an anharmonicity in the photonic spectrum: 

WQi^a /(a^i) ■ (6) 

Although not much is known for a realistic change of the photonic spectrum due to atoms, a possible 
anharmonic term is the interaction of the Hubbard model 

f{a''a) =uoa^a + U{a''af . (7) 

This quadratic anharmonicity probably overestimates the contribution of the atoms in the cavity. A more 
realistic case might be a sin^ anharmonic term: 

f{a)a) — ujQa)a + U sin^ (aa'''a) . (8) 

These two examples present a basis for analyzing the generic evolutionary properties of coupled anhar- 
monic microwave cavities. 



IV. TWO COUPLED CAVITIES 



The Hamiltonian of two coupled cavities 

H = -J{a\a2 + 4ai) + f{a\ai) + /(ajaa) (9) 

describes the tunneling of photons through the optical fiber with rate J. The properties of the two cavities 
are assumed to be equal and given by /(ajaj). While Fock states are eigenstates of the disconnected 

cavities (i.e. for Hq = f{a\ai) + /(a2a2)), this is not the case for the coupled system. This implies that 
a Fock state as the initial state undergoes an evolution by visiting the eigenstates of H . The average 
energy E is conserved, since the evolution of an isolated system is a unitary transformation: 

E = {-^tm^t) = (*o|-ff|*o) = fiN -k) + ,f(k) , (10) 

where k and N ^ k are the number of photons in the cavities for the initial state. 

The coupled microwave cavities with the two Fock states | TV, 0), |0, N) as initial states can be treated 
within the RPM. This is a generalization of the RPM of Ref. [10] for a single Fock state \N,Q) to a 
two-dimensional Fock basis. Assuming that N is even, all projected spaces 'H2j are two dimensional 
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and spanned by {\N — \j,N — j)} [j = 0, ...,N/2). This leads to a recurrence relation which can 
accommodate the Fock states |A^, 0), |0, N). In terms of Pauli matrices this reads {a'j}j=o,...,3 

gfc+i = Ofc+iCTo + 5fc+icri, go = ^ _ 2/( jV/2) ^° (A = 0, 1, iV/2 - 1) (11) 

with coefficients 

z - fk+i - J^ak{N/2 + k + l){N/2 - k) 



Ofc+l — 2 

z - fk+1 - .Pak{N/2 + k + l){N/2 - k) - J^hl{N/2 + k + lf{N/2 - kf 



jHk{N/2 + k + l)iN/2-k) 

^k+l - 72 



z - fk+1 - •Pak{N/2 + k + l){N/2 - k) - Pbl{N/2 + k + l)^N/2 - k^ 



(12) 



(13) 



and 



fk+i = f{N/2 + k + l) + f{N/2 -k-l). 
The initial conditions for the coefficients are 

, 1 

ao = Do 



z - 2f{N/2) 
and the termination of the iteration gives 

gN/2 = aN/2<^0 + bN/2(^l ■ (14) 

Here the Pauli-matrix structure is based on the two Fock states |0,-/V) and |-/V, 0). It should be noticed 
that the cavity mode ujq appears as ojqN in f^. Thus, this term leads only to a global shift of the energy 

Z — > Z — UJqN . 

From the iteration of the recurrence relations ()12|) . (fT3|) we obtain the coefficients aM/2 and 5jv/2 of 
gN/2- This allows us to evaluate the polynomials Pn{z), Qn+i{z) of the resolvent in Eq. ([2]) for the two 
Fock states \N, 0), |0, N) as well as for the entangled even and antisymmetric NOON states 

|7V007V) = -i=(|iV,0) + |0,7V)), \NOONy ^^{\N,0)-\0,N)) (15) 

as initial states. Eq. (fH)) implies for a Fock state as initial state 

{N,0\{z~H)-'\N,0) = {0,N\{z-H)-'\0,N)=aN/2 , (16) 
for the symmetric NOON state as initial state 

(iV00iV|(z - H)-^\NOON) = aN/2 + bN/2 , (17) 
and for the antisymmetric NOON state as initial state 

'(A^OOA^Kz - Hy^lNOON)' = aN/2 - bN/2 ■ (18) 

V. RESULTS 

We start with tunneling between two harmonic cavities, where f{N) = ujqN. In this case the denomi- 
nator of a7v/2> bN/2 factorizes: 

N/2 

QN+i{z)^zl[[{2kJ)^-z^] . (19) 

fe=i 
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FIG. 2: Photonic spectral function for 1000 photons in harmonic cavities with J = 1, e — 0.01, using symmetric 
and antisymmetric NOON states. The states with even parity are marked by crosses. 



In this expression we have dropped the global shift of the energy z by the harmonic term loqN. The 
distribution of the peaks of the PSD has a Gaussian envelope for both, the Fock and the NOON state as 
initial state (of. Fig. [2]). This can be directly explained by calculating the overlap between the eigenstates 
of the tunneling term |A^, k)t {k — 0, .., N) and the Fock state: 

(0,iV|iV,fc)t = 2-^/2(^ J . (20) 

Thus there is a binomial distribution for the overlap \ {0, N\N, k)t\'^ which becomes a Gaussian distribution 
for large N. The Gaussian result indicates a kind of Central Limit Theorem for independent photons. 
Such a behavior was also found previously for the free expansion of bosons from an initial Fock states 

For the case of anharmonic cavities the analytic expressions of the polynomials Pn{z), Qf^^i{z) are 
quite lengthy, and it is unknown if they can be written in a simple form such as Eq. ([19]). Therefore, it 
is more convenient to represent the results of the iterations for 1000 photons as plots of the PSD, using 
the expression of Eq. ([3|). The contribution of the harmonic term hujoN appears only as a global energy 
shift and is ignored in the subsequent discussion. 

The PSD is shown for harmonic cavities (Fig. [2]) as well as for anharmonic cavities with a Hubbard 
anharmonicity (Figs. [31 \^ and with a sin^ anharmonicity in Figs. [Sj [Sj The initial states are the two 
NOON states defined in Eq. ([T5|) . This allows us also to determine the PDS for the Fock states as initial 
states by linear combination. A characteristic feature of all the plots is the Gaussian envelope of the 
PDS, centered at the energy of the system in the absence of tunneling E — f{N). The tunneling rate J 
determines the width (cf. Fig. [4]). Moreover, the spikes of the PSD depend on the initial state, where 
the symmetric and the antisymmetric NOON state misses every other spike. This is also the case for 
harmonic and for anharmonic cavities. The reason is that the eigenstates appear either with even or with 
odd parity. Since the symmetric NOON state has even parity, its overlap with eigenstates of odd parity 
vanishes, and vice versa. The plots of the PSD in Figs. [3][ni indicate that every other state has even (odd) 
parity for all considered anharmonicities . The fact that states with odd and even parity are alternating 
results in a suppression of level crossing. 
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FIG. 3: Photonic spectral function for 1000 photons in Hubbard cavities with J = 1 and U — 0.001. The energy 
of the system is -E = 1000. 
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FIG. 4: Photonic spectral function for 1000 photons in Hubbard cavities with weaker tunneling rate J — 0.2 and 
U — O.QQl. The PSD is much narrower than for larger tunneling rates. 



VI. DISCUSSION AND CONCLUSIONS 



The overlap distribution function of the eigenfunctions with the initial state pj = \ {'$o\Ej)\'^ has been 
used in the literature to study various physical quantities, such as the entropy or the time evolution of 
the return probability of finite isolated quantum systems, e.g., the Bose-Hubbard model. This quantity 
is also useful to characterize the thermalization [16, 17j. For the Bose-Hubbard model this distribution 
function has a rich structure. For some parameter values it shows an exponentially decaying envelope 
though. In our plots the distribution function pj — K^'oliJj)^ is the envelope of the PSD, multiplied by 
e (cf. Eq. ©). 

Contrary to the exponential envelope of the Bose-Hubbard model in Refs. [la, we find a Gaussian 
envelope for the coupled cavities. For smaller systems we saw deviations from this behavior, at least for 
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FIG. 5: Photonic spectral function for 1000 photons in sin"^ cavities with U = W and J — 4, e — 0.01, a — 3.7119. 
The energy of the system is iS « 10. 
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FIG. 6: Photonic spectral function for 1000 photons in sin^ cavities with U 
same as in the previous figure. The energy of the system is iJ ~ 40. 



40, the other parameters are the 



the Hubbard anhamonicity [lO|. This structure is not very sensitive to the model parameter, as shown in 
Figs. BEl indicating a random-matrix type of behavior, similar to that found in other few-body systems 
[18l - [2(j | . In this context it is interesting to notice that a mean- field (Hartree-Fock) approximation of our 
model in Eq. ([9]) for a Hubbard anharmonicity gives a relatively simple behavior 15]: The photonic 
occupation Nt of one cavity is a periodic function of time t and reads 

Nt = ^[l + cn{Jt\NyN!)] , 

where the critical number of atoms is Nc — J/U. The Jacobian elliptic functions cn{Jt\N'^ /N^) is periodic 
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in the first argument Jt with equidistant frequencies 



En = vn, 
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f-iV2 sin2(e)/A^2 



1 



dQ . 



In particular, at the critical point N = Nc it decays like sech( Jt) = I/cosh(Ji). Such a behavior is quite 
different from our results of the Hubbard anharmonicity in Figs. [3l |4l 

The PSD reveals details of the evolution of photons which are tunneling between two microwave cavities. 
The fact that eigenstates appear pairwise with even and odd parity implies a strong dependence on the 
parity of the initial state. Thus the PSD provides a signature that the evolution of the return probablity 
|(\l/o|\E't)P depends on the choice of the initial state. This can be used to filter states according to parity 
in the evolution by choosing an appropriate initial state. 

The PSD is invariant with respect to the harmonic frequency of the cavities, except for a global 
energy shift. This reflects an important universality of the PSD that allows us to separate the harmonic 
from the anharmonic properties of the cavities. The behavior of the PSD for harmonic and anharmonic 
cavities is characterized by a Gaussian envelope, implying that the overlap of the eigenfunctions with 
the initial state pj = |(\E'o|£'j)P has a discrete Gaussian distribution. This simple feature is valid for 
sufficiently large photon numbers with deviations for smaller photon numbers though. Comparing the 
results of the RPM with earlier results of a Hartree-Fock approximation indicates that the latter is not 
a good approximation of the PSD, since the energy levels are very different. This does not rule out that 
Hartree-Fock approximation can be a useful tool for expectation values of some physical quantities. 
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